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equations exphcitly solvable in terms of orthogonal polynomials or associated special 
functions. 
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1. Introduction 

Many problems in quantum mechanics and mathematical physics lead to equations of 
the type 

a{s)y"{s) + T{s)y'{s) + \y{s) = (1) 

where a{s) and r(s) are polynomials of at most second and first degree, respectively, 
and A is a constant. These equations are usually called equations of hypergeometric type 
[To], and each of them can be reduced to the self-adjoint form 

[a{s)g{s)y'{s)]' + Xg{s)y{s) = (2) 

by choosing a function g such that (o'q)' = tq. 

The equation ([T]) is usually considered on an interval {a,b), chosen such that 

(j(s) > for all s G (a, b) 

g{s) > for all s G (a, b) (3) 
\ims^a(j{s)g{s) = lims^b (J is)g{s) = 0. 

Since the form of the equation ([1]) is invariant under a change of variable s ^ cs + d, 
it is sufficient to analyse the cases presented in table 1. Some restrictions are to be 
imposed to a, (3 in order the interval (a, b) to exist. 



Table 1. The main cases 





r(s) 


q{s) 


a, 13 


(a, 6) 
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as + P 




a < 
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as+[3 




a < 0, /3 > 


(0,oo) 


1-^2 


as + P 


(l + 5)-(a-/3)/2-l(^l __5)-(a+/3)/2-l 


a < 13 < —a 


(-1,1) 




as + P 


(s + l)("-/3)/2-l(5_l)(a+/3)/2-l 


-/3 < a < 


(l,oo) 


s^ 


as + P 




a < 0, /3 > 


(0,oo) 


s^ + l 


as+(3 


^■j^ _j_ ^2^a/2— lg/3 arctan s 


a < 


M 



The equation ([T]) defines an infinite sequence of orthogonal polynomials in the case 
cr(s) G {1, s, 1 — s^}, and a finite one in the case 0"(s) G {s^ — 1, s^, + 1}. The 
literature discussing special function theory and its application to mathematical and 
theoretical physics is vast, and there are a multitude of different conventions concerning 
the definition of functions. The table 1 allows one to pass in each case from our 
parameters a, /3 to the parameters used in different approach. 

In section 2 we briefly present some results concerning orthogonal polynomials, 
associated special functions, hypergeometric type operators and related Schrodinger 
type operators which are needed in section 3. In quantum mechanics there exist 
potentials, called quasi-exactly solvable, for which it is possible to flnd a flnite portion 
of the energy spectrum and the associated eigenfunctions exactly and in closed form 
[D El m [HI [m [12]. An algorithm for generating new explicitly solvable systems and 
some applications are presented in section 3. 
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2. Orthogonal polynomials and associated special functions 



Let t{s) = as + P he a fixed polynomial, and let 



A, 



r 



'{s)i- 



(4) 



2 ' ' ' ' 2 
for any £ G N. It is well-known that for A = , the equation ([T]) admits a polynomial 
solution = of at most i degree 



0. 



If the degree of the polynomial $£ is i then it satisfies the Rodrigues formula [TU] 



B, d^ 



[a\s)Q{s)] 



(5) 



(6) 



(7) 



A: 



(8) 



q{s) As^ 

where is a constant. Based on the relation 

{ (5 G M I lims_5.a o"(s)f)(s)s'^ = lims^;, cr(s)f)(s)s^ = } 

[0,cx)) if a{s) G {1, s, 1 - s^} 

[0,-a) if a{s) G {s^ - 1, s^, + i} 

one can prove [21 |3l [10] that the system of polynomials {$£ | £ < A}, where 

oo for a{s) G {1, s, 1 - s^} 

ifi for a{s) G {s^-l, s^, s^ + l} 

is orthogonal with weight function q{s) in {a,b). This means that equation ([T]) defines 
an infinite sequence of orthogonal polynomials 

*0, $1, '^'2, ••• 

in the case a{s) G {1, s, 1 — s^}, and a finite one 

$0, $1, $L 

with L = max{£ G N | £ < (1 - a)/2} in the case a{s) G {s^ - 1, s^, + 1}. 

The polynomials can be expressed (up to a multiplicative constant) in terms 

of the classical orthogonal polynomials as [3] 

Hei^J^s-^ 



2 

as J 



-2a 



p{(a-/3)/2-l, (a+/3)/2-l),_^^ 



l-a-2i 



(!) 



in 


the 


case 


a 




= 1 




in 


the 


case 


a 




= s 




in 


the 


case 






= 1 - 


s' 


in 


the 


case 






= ^2- 


- 1 


in 


the 


case 


a 




= S^ 




in 


the 


case 


a 


's) 


= s'^ 


- 1 



(9) 



jfp«"+i/3)/2-l, (a-i/3)/2-l)^..^^| 

where Hg, and P^^'''^ are the Hermite, Laguerre and Jacobi polynomials, respectively. 
The relation iQ does not have a very simple form. In certain cases we have to consider 
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the classical polynomials outside the interval where they are orthogonal or for complex 
values of parameters. 

Let ieN, i<A, and let me {0, 1, £}. The functions 



^U\s) = «^"^(^)^*r^(^) Where = (10) 

are called the associated special functions. If we differentiate ([5]) m times and then 
multiply the obtained relation by k"^{s) then we get the equation 

n^^^;f^ = A,$(-^) (11) 

where Hm is the differential operator 

^r. = -ais)i-ris)i + ^^^ 

+ - l^i^ - 2)^"(s) - mr'is). 



7{S) 

For each m<A, the special functions with m<i<A are orthogonal with respect 

to the scalar product 



{f,9)= / f{s)g{s)g{s)ds. (13) 

J a 

The operators Hm are directly related to some Schrodinger type operators. If 
(a, b) — >■ (a', b') : s X = x{s) is a differentiable bijective mapping such that 
dx/ds = ±1/k(s) and {a',b') — )■ {a,b) : x t-j- s{x) is its inverse then the functions 

= V<s{x))g{s{x))^^;:;^'\s{x)). (14) 

with m < i < A are orthogonal [2] 

£'^t'^(^)^i'^'H^)rf^ = for ^^fc 
and satisfy the equation [21 S] 
rf2 



+ Kn(x) 



<'^ = A.<'^ (15) 



where Vm{x), defined in terms of the function 77(5) = l/yK(s) g{s), is given by 
VUx) = + - |m(m - 2)a"(.) 



(16) 



J s=s(x) 

For example, in the case a{s) = 1, the change of variable M— j-M : x 1— )■ s{x) =x leads to 
V^x) = ^x' + ^x + ^ + ^-am. (17) 
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3. New explicitly solvable Schrodinger type equations 

In the case of a second order differential equation 

^(r)|^ + 5(r)|: + C(r)]^(r) = 
with A{r) 7^ we can eliminate the first order derivative by using the function 

h{r) = exp (^f ^^dt^ . 
The equation ffTSl) is equivalent to the equation 

^ h{r) [A{r)^ + B{r)i + C{r)\ ^ /.(r)^(r) = 
which can be written as [5l [9] 



(P _|_ AA{r) C{r)-2A(r) B'(r)+2B{r) A'{r)--B^{r) 



The Schrodinger type equations ( IT5|) have the form [1] 



h{r)tlj{r) = 0. 



(a;) = 0. 



(19) 
(20) 
(21) 
(22) 



If, for k^{—l, 1}, there exists a differentiable bijective mapping 
(a, 6) — )■ (a', ft' ) : r i-)- x(r) 

such that 

x'{r) = 



{a',b') 
1 

then the equation (l22l) is equivalent to 

■^(^_kK^,p,in) /fe(x(r)) 16 (/fc(x(r)))3 



-E 



A/4RO)*5"i^(^W)=o 



(23) 



' 4 (/fe(a;(r)))^ 

where = — Cfc(a, /3, m). We know that a, /3, £ and m must satisfy certain restrictions. 
In |5] the authors consider the Schrodinger equation (translated harmonic oscillator) 



-^+e''x' + px + x 



(x)=0. 



In this case the substitution x = \/1r leads to the equation 



(P- I P \ A ^ J_ _L 



3 1 



dr'^ ' V2r ' 2r- 16 

and the substitution x = ^ (3r/2)2 to the equation 



^0(v/2f)=O. 



' dr^ 



_5__1 

36 r 



3+P 



?0(^(3r/2)^ 



:0. 



(24) 



(25) 



(26) 



In [5] these equations are considered as two new exactly solvable Schrodinger equations. 
This is not obvious because, in equation flMj) . the eigenvalue A depends on 6*^ and p. 
The Schrodinger type equation corresponding to the potential (E 



dx^ 



-^x^ + ^x + ^ + %-am + al 



(x) = 



(27) 
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is satisfied for a < and any m, £ G Z sucli tliat 0<m<i. It is similar to 

I^i{x)=x, Ci{a, f3,m) = C-i{a, f3,m) = ^ and Co{a, f3,m,i) 

Using tlie substitution x = ^ (3r/2)2 = 



hix) = x'^, 



rs we get tlie equation 



and by using tlie substitution x = y2r the equation 



1 3_ J_ 

2r 16 r2 



2r)=0. 



Let Ci and C2 be two fixed real numbers. The function 

iPe,^ : (0, oo) M, V^,,„(r) = ^ ^^"iH^rf; 
is a solution of the equation 



dr2 +"^1 ( ^2" ) + '^2 ( ) ^ 36 ^2 



il) = Eijj 



for = — ^ if a, /3, £, m satisfy the system 



Cl 



^ + I - am + a£ = C2 
and the conditions m,iE{0, 1, 2, ...}, m<i, a<0. In the case ci >0, the functions 



satisfy the relation 



{-2v/cr , ±2^C2+v^(l+2€-2m)) 



* 



± 



± 



for 



= ±2 ^/ci C2 + Cl 751(1 + 2^ - 2m 
if m, £g {0, 1, 2, ...} are such that 
m < £ 



(28) 

(29) 
(30) 
(31) 

(32) 

(33) 
(34) 
(35) 



and C2 + ycr(l + 2£ - 2m) > 0. (36) 
By using ([9]), ( |TOl) . (fT4l) and table 1, the solutions ipf^ can be expressed in terms of 



Hermite polynomials 



re 



r exp 



3 3/3 



y^ra ± + v^(l + 2e- 2m) r 



X 



T ^ Vc2 + v^(l + 2£- 2m] 



(37) 



and one can remark that they are square integrable on (0, 00). In view of the well-known 
relation H'^ = 2i H£_i, the function ^f^ir) coincides up to a multiplicative constant to 
the function 

^±(r) = ^ exp (^-lyiy^rl ±0^C2 + ^{1 + 2n) 
X Hn (^x T ^ + v^(l + 2ra) 



2 

r 3 



(38) 



r3 
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where n = i — m. For Ci >0, the function ip^i''^) eigenfunction of the Schrodinger 

type operator 

^ = -S^ + c,{fr + c,{i,y-U (39) 
corresponding to the eigenvalue 



= ±2y/ci C2 + civ^(l + 2n) (40) 
for any n G {0, 1, 2, ...} satisfying the relation 

n> (41) 

Certain solutions of the equation (see (!29|) ) 

^ = Ei: (42) 
may be found by looking for solutions of the system of equations 



a/3 

q2 



2 (43) 



^ + I - am + a£ = C2 

satisfying the conditions a <0 and m, £g {0, 1, 2, ...} with < m < i. The corresponding 
value of is — The solution of the system leads to an equation of third degree, 
the formulas are more complicated and will not be presented here. The other cases 
presented in table 1 may also lead to some new explicitly solvable systems. 



4. Concluding remarks 

The method proposed by Derezihski and Wrochna allows us to generate new explicitly 
solvable systems if we take into consideration the parameter dependence of the 
eigenvalues of the starting system. Some interesting explicitly solvable systems can 
be obtained by starting from solvable hypergeometric type equations containing several 
parameters. 
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